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Meson properties at finite temperature and density are studied in lattice QCD simulations with 
two-flavor Wilson fermions. For this purpose, we investigate screening masses of mesons in 
pseudo-scalar (PS) and vector (V) channels. The simulations are performed on 16^ x 4 lattice 
along the lines of constant physics at mps/m^/\T=o = 0.65 and 0.80, where mps/mv|r=o is a ratio 
of meson masses in PS and V channels at T = 0. A temperature range is T /T^c = (0.8 —4.0), 
where Jpc is the pseudo-critical temperature. We find that the temperature dependence of the 
screening masses normalized by temperature, Mq/T, shows notable structure around Tpc, and 
approach 27t at high temperature in both channels, which is consistent with twice the thermal mass 
of a free quark in high temperature limit. The screening masses at low density are also investigated 
by using the Taylor expansion method with respect to the quark chemical potential. We find that 
the expansion coefficients in the leading order become positive in the temperature range, and 
thermal and density effect on the meson screening-masses becomes apparent in the quark-gluon 
plasma phase. The meson screening-masses are also compared with the gluon (Debye) screening 
masses at finite temperature and density. 
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1. Introduction 

Study of hadron properties at finite temperature and density is important to understand behav- 
ior of quarks and gluons in hot and/or dense QCD medium. In particular, mesons in the medium are 
expected to have abundant information about characteristic properties of QCD, such as the decon- 
finement phase transition at finite temperature and the partial chiral symmetry restoration at finite 
density In this article, we focus on screening masses of mesons in pseudo-scalar (PS) and 
vector (V) channels calculated from the spatial correlation functions, and present current results 
of our group in lattice QCD simulations with two flavors of the Wilson-type quark action. Us- 
ing gauge configurations generated by WHOT-QCD Collaboration [^, ^], we calculate the meson 
correlation functions along the line of constant physics and extract temperature dependence of the 
meson screening-masses from long spatial-distance behavior of the correlators. 

We find that the screening masses in both channels show notable behavior around pseudo- 
critical temperature Tpc, where the transition occurs from the hadronic phase to the quark-gluon 
plasma phase. On the other hand, these approach InT at high temperature, which corresponds to 
twice the thermal mass of a non-interacting quark in high temperature limit. We also investigate 
the meson screening-masses at low density by using the Taylor expansion method with respect to 
the quark chemical potential. We find that the leading order expansion coefficients become positive 
in the temperature range, and thermal and density effect on the meson screening-masses becomes 
apparent in the quark-gluon plasma phase. We also compare the meson screening-masses with 
the gluon (Debye) screening-masses obtained from Polyakov-line correlation functions, and find 
characteristic difference in medium contributions between mesons and gluons. 

The paper is organized as follows. In Sec. 2, we show the formalism to calculate the screening 
masses on the lattice. In Sec. 3, numerical results of the screening masses at finite temperature 
and zero density are discussed, and results of the screening masses at finite density in Sec. 4 using 
the Taylor expansion method. In Sec. 5, we compare the meson screening-masses with the gluon 
(Debye) screening masses. Section 6 is devoted to summarize the paper. 

2. Meson screening masses on lattice 

In order to extract the meson screening masses at finite temperature and density, we calculate 
an expectation value of the spatial correlation functions G{x) of the mesons, 

{G{x)) = £(M(x,>',z,0M(0, 0,0,0)^), (2.1) 

where M{x,y,z,t) = q{x,y,z,t)rq{x,y,z,t) is the meson operator with F denoting the gamma ma- 
trix, i.e., r = 75 for a pseudo-scalar (PS) meson and F = 7^ for a vector (V) meson. The correlator 
is summed over y,z,t, which means the zero-momentum projection for y and z directions and zero- 
energy projection for temporal direction. 

In order to extract response to finite density, we apply the Taylor expansion method with 
respect to the quark chemical potential, n = + jUrf)/2, where /i„ (/i^/) is the chemical potential 
for the u (d) quark. The Taylor expansion enables us to investigate meson properties at /i 7^ 
from the expectation values at /i = 0. Following Ref. we expand the expectation value of the 
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operator in powers of /i = /i/T as, 

/ &U f'"'^!!'"™(detD(/x))2 



/ &U f'-'^stan(detD(/x))2 



J m e-^^'"-" (A + A/t + 5 A^t 2 + 0{p.^)) 

+ + i^A' + o(A^))(i + |a + + ojfi'))) 
+ ((^^> + A + ((^^) + \m + ^(^^) - + ^(^')' ^2-2) 

where A = (detD(/x))^ with denoting the Dirac operator, and the dots on operators denote the 
derivatives with respect to p.. Note that the expectation values in the right-hand-side of Eq. ( p!2| ) 
are calculated at = 0. We take to be the meson correlator G = tr(D^' (jU„)rDQ^^ (;U^)r^) where 
tr denotes the trace with respect to color, spinor, and flavor indices. Then the Taylor expansion of 
the meson correlator is given by, 

(G)^ = (G)o + (G)iAt + (G)2A' + O(A'), (2.3) 
where the expansion coefficients become, 

(G)o = (tr[D;oir75(D-i)|o75r^]), (2.4) 
(G)i = 0, (2.5) 

(G)2 = (Gopr)2+(Gdet)2, (2.6) 
(Gopr)2 = + 

= 2(Retr[(D-iDD-iDD-i),or75(D-i)^75rt]) 

-(Retr[(D-iDD-i)^r75(D-i)]o75r^]) 

-(Retr[(D-'Z)D-i).ri,r75(D-iZ)D-i)^75r^]) 

+4(Imtr[(D-iDD-i),or75(D-i)]o75r1^].ImTr(D-iD)), 

1 , A, 1 , , ,A, 
(Gdet)2-2^G-)--(G)(-) 

= Re{(tr[D;oir75(D-i)Jo75r1-](2(Tr(D-iD))2 -Tr(D-iDD-'D) +Tr(D-iD))) 
-(tr[D;oir75(D-i)Jo75rt])(2(Tr(D-iD))2-Tr(D-^DD-^D)+Tr(D-'D))}, 

where Tr denotes the trace including space-time indices in addition to those for tr. We have divided 
the second derivatives into two parts: (Gopr)2 and (Gdet)2> where the former includes the derivatives 
of the operator G, and the latter consists only of the derivatives of the quark determinant. Note that 
the meson correlator does not have the odd orders in the Taylor expansion since it is symmetric 
under jj. — )• — /i, namely the meson correlator is invariant under the charge conjugation. For the 
calculation of the trace with respect to spatial indices, we apply the random noise method with 100 
sets of U(l) random numbers. 
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In order to study the screening effect, we fit tfie meson correlator by the following form, 



(G(x))(^,r)=A(At,r) 



'^-M{nJ).x _^ ^-M(Ai,r)(L--v)l ^ (-2.7) 

where L is the spatial lattice size, and we assume that contributions of finite jJ. appear only in the 
coupling factor (A) and the meson screening mass (M). We also assume that A(/i, T) and M{ix,T) 
are also expressed as power series in fx, 

A{^,T) = Ao +A2fL^ + O(A^), (2.8) 

M(At,r) = Mo+M2A^ + 0(/i^). (2.9) 

By comparing both sides of Eq. (2.8) at each order of /i, we obtain, 

(G(x))2 



A^Ae-'^^' + e-^'^'^^-'^' 



^ 1 tanh 



(2.10) 
(2.11) 



{G{x))o 

We extract the screening masses by fitting the meson correlators for each temperature with the 



expressions Eqs. (2.1C) and (2.11) at large distance. 



3. Numerical simulations 

Simulation setup is the following. We utilize the gauge configurations generated by WHOT- 
QCD Collaboration on 16^ x 4 lattice with the renormalization-group improved Iwasaki gauge ac- 
tion and Nf = 2 clover-improved Wilson quark action The simulations have been performed 
along the line of constant physics corresponding to the PS and V meson mass ratio, mps/mv|r=o = 
0.65 and 0.80 at T = 0. The temperature range for mps/mv|r=o = 0.65 (0.80) is T /Tp^ =0.82^.0 
(0.76-3.0), where Tpc is the pseudo-critical temperature for the transition from hadronic phase to 
quark-gluon plasma phase. The number of configurations we use is 100 for each temperature and 
quark mass. 

3.1 Screening masses at finite temperature and zero density 

Figure 1(a) shows the meson screening-masses normahzed by temperature, Mq/T, in PS chan- 
nel as a function of T /Tpc- The circle (triangle) points correspond to the results at nips/mv|r=o = 
0.65 (0.80). The same figure in V channel is shown in Fig. 1(b). We can see a concave structure 
around Tpc, that is, when temperature increases, Mq/T decreases below Tpc, whereas it increases 
above Tpc. This implies that the screening masses Mq stay constant below Tpc, whereas they mono- 
tonically increase above Tpc. Namely, thermal effect on the meson screening-masses becomes 
apparent in the quark-gluon plasma phase. At high temperature, Mq/T converges to a constant 
value of 271, which implies that the meson becomes a weakly interacting pair of a quark and an 
anti quark, each carrying the thermal mass nT. We also find clear quark-mass dependence, i.e., 
magnitude of the screening masses with lighter quark mass (mps/mv|r=o = 0.65) becomes smaller 
than that with heavier quark mass (mps/mv|r=o = 0.80) in both channels, similarly to the ordinary 
meson mass measured by temporal correlation. This implies that the screening of the meson spatial 
correlation becomes weak when the quark mass becomes small. 

The temperature dependence of the meson screening masses is consistent with that calculated 
in the staggered-type quark action 
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Figure 1: Meson screening masses Mq/T in the PS channel (a) and V channel (b) as a function of tempera- 
ture. 



3.2 Screening masses at finite density 

In this section, we investigate properties of the screening masses at finite density by calculating 
the second response of the screening masses M2/T to the quark chemical potential jJ. via the Taylor 
expansion method. Figure 2(a) and 2(b) show temperature dependence of M2/T in PS and V 
channels, respectively, with mps/mv|r=o = 0.65 (circle points) and 0.80 (triangle points). We 
find that M2/T is always positive in the temperature range we have explored, which implies that 
the screening masses increase in the leading order contribution of jx. M2/T increases rapidly at 
Tpc which means that density effect on the meson screening-masses becomes significant in the 
quark-gluon plasma phase. We also found from the simulations that main contribution to the rapid 



increases of M2/T comes from the (Gopr)2 in Eq. ( |2.6| ) which consists of the derivatives of the 
operator. We can also see the quark-mass dependence in the PS channel; the magnitude of M2/T 
with lighter quark mass (mps//Mv|r=o = 0.65) is slightly smaller than that with heavier quark mass 
(mps/mv|r=o = 0.80). This is similar to the quark-mass dependence of Mq/T, but different from 
the results of M2 calculated in the staggered-type quark action [Q]. This difference should be further 
investigated in the simulations with smaller quark mass and larger spacial volume. 



4. Comparison with gluon screening-masses 

Let us compare the screening masses of the PS mesons with these of the gluons. The response 
to the quark chemical potential jj. should be different for fermionic objects and gluonic ones, be- 
cause the former (such as the meson screening-mass) have direct coupling with /i described by 
(Gopr)2 while the latter (such as the gluon screening-mass) have only indirect coupling via the dy- 
namical quark loops in the medium. Therefore, it is important to compare these two screening 
masses. The gluon screening-mass (so called Debye screening mass) has been studied from the 
Polyakov-line correlator based on the Taylor expansion method using the same gauge configura- 
tions [|,|]. The Debye screening mass is also expressed as power series of p., 

mD{}X)=mDfl + mD,2Pi^ + 0{}jL^). (4.1) 
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Figure 2: Second response of the screening masses M2/T with respect to /i as a function of T/Tpc for PS 
channel (a) and V channel (b). 
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Figure 3: (a) meson screening-masses in PS channel, (b) gluon (Debye) screening-masses at mps/mv\T=o - 
0.65 as a function of T/Tpc- 



Figure 3 shows the meson (a) and Debye (b) screening-masses at mps/mv\T=o = 0.65, where 
the circle (triangle) plots show Mq/T and mjj o/T {M2/T and niD^/T), respectively. We find 
characteristic difference in the temperature dependence, i.e. the meson (Debye) screening-masses 
increase (decrease) when temperature increases above T^c- This is related to the fact that, at high 
temperature limit, Mq goes to 2nT which is twice the thermal mass of a free quark, whereas wio 
is proportional to the running coupling and goes to zero according to the prediction of the thermal 
perturbation theory. We also find that the ratio of the meson screening-masses, M2/M0, is larger 
than that of the Debye screening-masses, 2/"^d,o- This means that the response of the fermionic 
object to n is larger than that of the gluonic object due to the difference in their ways of coupling 
with jj. mentioned above. 

5. Summary 

We have studied the meson screening-masses in PS and V channel at finite temperature and 
density in lattice QCD simulations with two-flavor Wilson fermions. The simulations have been 
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performed along the line of constant physics at mps/mv\T=o = 0.65 and 0.80 with the temperature 
range of T /Tpc = 0.82^.0 and 0.76-3.0, respectively. On the basis of the Taylor expansion method, 
we have calculated temperature dependence of the leading order term of the screening masses (Mq) 
and the second response to the quark chemical potential /i (M2). We have found that Mq/T shows 
a concave structure around Tpc and goes to 271 at high temperature, which corresponds to twice the 
thermal mass of a free quark. From the quark-mass dependence, we have seen that the screening 
effect on the meson correlator becomes weak when the quark mass becomes small. 

The second response M2/T is always positive in the temperature range we have explored, 
which implies that the screening masses increase in the leading order contribution of jj.. We have 
also found that thermal and density effect on the meson screening-masses becomes significant in 
the quark-gluon plasma phase. 

We have also compared the meson screening-masses with the gluon (Debye) screening-masses, 
and found characteristic difference of the temperature dependence, i.e. the meson (Debye) screening- 
masses increase (decrease) when temperature increases. We have also seen that the response to jJ. is 
different between the meson and Debye screening-masses, which reflects the fact that the fermionic 
object has direct coupling with /i, whereas the gluonic object has only indirect coupling via the dy- 
namical quark loops in the medium. 
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